subroutine r2abij_t2abmn_wmnij(r2, t2, &
              wmnij, icore, fact)
!
use mod_ioff
use mod_orbit 
implicit none
!
real*8, intent(in) :: fact
real*4, intent(in) :: wmnij(*)
real*4, intent(in) :: r2(*)
!
real*4, intent(inout) :: t2(*), icore(*)
!
integer :: i, j, a, b, nsize
real*8, external :: dnrm2
!this is to do r2 = r2 + 0.5*t(ab,mn)*Wmnij
! AAAA 
call VpqmnVmnrs_to_Vpqrs(t2(ioi3(1)), wmnij(ioi1(1)), 1, 1, popa, popa, 0, &
               vrta, vrta, 0, popa, popa, 0, r2(ioi3(1)), 1, fact*2.d0, 1.d0)
! ABAB  
call VpqmnVmnrs_to_Vpqrs(t2(ioi3(2)), wmnij(ioi1(2)), 1, 1, popa, popb, 1, &
               vrta, vrtb, 1, popa, popb, 1, r2(ioi3(2)), 1, fact*2.d0, 1.d0)
! BBBB  
call VpqmnVmnrs_to_Vpqrs(t2(ioi3(3)), wmnij(ioi1(3)), 1, 1, popb, popb, 0, &
               vrtb, vrtb, 0, popb, popb, 0, r2(ioi3(3)), 1, fact*2.d0, 1.d0)
!
return
end
